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Since the early observation of nematic phases of disc-like clay colloids by Langmuir in 1938, the 
phase behaviour of such systems has resisted theoretical understanding. The main reason is that 
, there is no satisfactory generalization for charged discs of the isotropic DLVO potential describing the 

effective interactions between a pair of spherical colloids in an electrolyte. In this contribution, we 
show how to construct such a pair potential, incorporating approximately both the non-linear effects 
of counter-ion condensation (charge renormalization) and the anisotropy of the charged platelets. 
' The consequences on the phase behaviour of Laponite dispersions (thin discs of 30 nm diameter and 

1 nm thickness) are discussed, and investigation into the mesostructure via Monte Carlo simulations 
are presented. 



I. INTRODUCTION 



o 
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With the possibility to form an orientational ordering, the phase behaviour of anisotropic colloids is richer than its 
counterpart for spherical particles. Surprisingly, the isotropic-nematic transition expected on purely entropic grounds 
(excluded volume effects) has been extensively reported for rod-like colloids in the last sixty years |ij , but could only 
be observed recently for discotic particles in apolar media This experimental work directly points to the subtle 
q ■ effect of electrostatic interactions, given that in the widely studied model system of colloidal platelets, namely aqueous 
l — l [ clay dispersions, the isotropic/nematic phase separation is hindered by a ubiquitous "fluid-solid" transition. 

In spite of an important experimental and theoretical effort in the last 10 years, triggered by the emergence 
of Laponite as a model system for disc-like colloidal suspensions ||, the above transition is far from being well 
understood PH3l|. In particular, the precise nature of the phases observed experimentally is debated ||,|l4|,g4| , so 
(---. that the relevance of the traditional terminology "sol-gel" to describe the transition is somewhat controversial. 

On the other hand, theoretical approaches to describe the phase behaviour of charged disc-like particles in an 
C ■ ' electrolyte are in their infancy [pp|, p7| , ^ | , essentially because there is no satisfactory generalization for discs of the 
isotropic DLVO potential |3^] describing the effective interactions between a pair of spherical colloids. In the simplest 
situation of two coaxial and parallel uniformly charged platelets, the effective Coulomb force has been computed within 



linear 
out |33[ 



ll| and non-linear [^l] Poisson-Boltzmann (PB) theory. A screened electrostatic pair potential has been worked 
allowing to compute analytically the interaction energy for discs of arbitrary orientation; this approach holds 
at the level of linearized PB theory (weak electrostatic coupling) , and is expressed as a perturbative expansion in the 
parameter nro, ro being the disc radius and k the inverse Debye length in the electrolyte. It consequently becomes 
less reliable when ktq becomes of order 1 or larger (the typical situation for Laponite solutions), or when the charge 
on the platelets becomes too large (high electrostatic coupling, also typical of clay systems). The latter limitation 
may be circumvented by introducing the concept of charge renormalization |34]-[56],[32| , while the former requires the 
re-summation of all the powers of ktq involved in the "multi-polar" expansion propounded in |3q |. In this work, we 
discuss how to take both aspects into account and investigate the consequences on the phase behaviour. Our goal is 
to understand the effect of colloid anisotropy on a possible phase transition. We show that considering electrostatic 
effects only and discarding van der Waals interactions yields qualitative differences between spherical and discotic 
colloids and that the equivalent of the fluid-solid transition for spheres is nevertheless reminiscent of the "sol-gel" 
transition reported in ||. 

The paper is organized as follows. In section ||, we obtain a pair potential for discs valid at large distances within 
linearized PB theory. We show that at asymptotically large separations, this potential remains strongly anisotropic, 



unlike its counterpart in vacuum or in a plain dielectric medium (no salt). The approach is generalized in section III 



to an arbitrary "one-" or "two-dimensional" colloid (i.e. with vanishing internal volume). The non- linear phenomenon 



of counter-ion "condensation" leading to a renormalization of the bare platelet charge is then considered in section IV 
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and the resulting effective charge successfully tested against "exact" numerical simulations reported in the literature. 
In section [v|, we propose a first and simplified investigation into the phase behaviour by mapping the interaction energy 
onto an isotropic Yukawa potential. Finally, the full anisotropic potential including charge renormalization is used 
to investigate the mesostructure of an assembly o f int eracting discs by means of Monte Carlo computer simulations 
(section VI), and conclusions are drawn in section VII. 



II. LARGE DISTANCE SCREENED ELECTROSTATIC POTENTIAL 



We consider first a unique platelet V of radius rg and surface charge a in a 1:1 electrolyte of bulk density n 
(infinitc dilution limit). The solvent is assumed to be a dielectric continuum of permittivity e. Within linearized 
Poisson-Boltzmann theory (LPB), the dimensionless electrostatic potential obeys the following Poisson's equation 

V^^K^, (1) 

where n 2 = %-ni B n is the inverse Debye length squared and Ib — e 2 /(kTs) denotes the Bjerrum length corresponding 
to the distance where the bare Coulomb potential felt by two elementary point charges becomes equal to the thermal 
energy kT (£ B ~ 7A in water at room temperature). The potential is chosen to vanish at infinity and e denotes the 
elementary charge. 

We are interested in the behaviour of <f> at large distances, which is relevant to describe the interactions in dilute 
suspensions. The solution of Eq. (|l]) has been worked out in the form of an integral representation in p^Jll]] , or as a 
multi-polar- like expansion fl33[ . It is however more convenient to obtain an explicit expression, that is a good starting 
point to derive the pair interaction. This can be achieved by writing the solution as a convolution of the surface 
charge of the platelets with the screened Coulomb potential: 

, / x „ f cr exp(— k|t — s|) , 9 

* r > = He V-b| ^ (2) 

It is important to realize that such a convolution procedure would give incorrect results for polyions of non vanishing 
excluded volume (e.g. spheres or rods of non zero radii, see the appendix). Introducing the unit vector n normal to 
the disc surface, and the angle 9 £ [0,7r/2] between the corresponding direction and the position vector r with origin 
at platelet centre, we obtain the leading contribution 

s «r»i „ Ii(/cr o sin0) e~ Kr 2 

4>{r) ~ £ B Z b aro2 — with Z barc e = nr a : (3) 

nro sin 9 r 

where T denotes the modified Bessel function of the first kind, such that \\{x) ~ x/2 for \x\ « 1. A peculiarity 
of screened Coulomb potential appears at this point: whereas the large r potential becomes isotropic and behaves 
as Zbaxe/r m vacuum (situation corresponding to the limit k — > 0), the anisotropy is present at all distances in an 
electrolyte (k ^ 0): the 9 and r dependence factorize in equation (0). In practice, the anisotropy of the potential is 
generically significant when the size of the object under consideration is larger than the Debye length (i.e. ktq > 1 
here) . For Laponite discs of radius 150 A this crossover corresponds to an ionic strength 7* = 10~ 4 M. As noted 
in |l8|| , /* corresponds experimentally to a threshold value delimiting qualitatively different phase behaviours. At 
fixed distance r, the potential (|3|) is minimum for 9 = which corresponds to the configuration where the "average" 
distance between the point where <f> is computed and the platelet is maximal; on the other hand, the potential is 
maximum for 9 = it/2. From (^ the ratio between these extremal values is r-independent, and reads 2Ii(Kro)/(«Xo) 
which can be as large as 10 for ktq = 5 (important anisotropy). 

This result may be used to compute the potential energy of interactions between two platelets A and B with 
arbitrary relative orientations, as shown in Fig. This energy is obtained by integrating the screened potential 
created by disc A defined by Eq. (||) over the surface charge distribution of the second disc (B): 

V AB (r,9 A ,9 B ) = [ -0 A (r + s)d 2 s, (4) 
Jb e 

where Vab is the dimensionless potential expressed in kT units. Again, such a procedure would not give the proper 
interactions in the case of polyions with internal volume. At large distances, we obtain 

V AB (r,0M ~ Ob* ^ 3 T^osinM e^: 

ktq sin 9 A Kro sin 9 B r 
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Clearly, the relative orientation of the platelets is not completely specified by the three parameters r, 9 a and Ob, but 
the omitted Euler angles only appear in higher order terms such as exp(— «x)/r 2 , exp(— nr)/r 3 . . .In particular the 
energy (||) appears insensitive to a precession of the discs around their centre to centre direction. 

At fixed centre to centre distance, the above energy (always repulsive |3q| ) is maximized for co-planar discs (9 a = 
Ob = t/2, which corresponds to the maximum overlap of electric double layers) and minimized when the discs are co- 
axial and parallel (Oa = Ob = 0, see the configuration represented in Fig. ||). A situation of intermediate electrostatic 
energy is that of T-shape perpendicular discs (0a = and Ob = t/2). These results are consistent with the numerical 
linearized Poisson Boltzmann pair potential reported in ^}\ , where it was also shown that the interactions between 
platelets at constant surface charge or constant potential were qualitatively very similar. 

In reference p3| , the potential energy Vab was obtained in the form of a multi-polar expansion. Unlike its unscreened 
counterpart where the multipolc of order I has a l/r l+1 large distance contribution to the potential, this expansion 
is such that the multipole of order I contributes to orders exp(— nr)/r, exp(— Kr)/r 2 . . . exp(— nr)/r l+1 . The large 
r potential can in principle be obtained by re-summation of all multi-polar contributions, which should lead to 
expression (|J). Alternatively, truncating the multi-polar expansion at a given order amounts to expanding (||) in 
powers of ktq. Including monopolc-monopole, monopole-quadrupole and quadrupole-quadrupole interactions, the 
following expression was obtained in [p3| : 



Vab(t, a ,( 



= Z, 



bare 



(sin 2 



4 4 
Kit 
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^■(sin 2 9 A sin 2 9 B ) 



(6) 



From (|5|), we obtain 
V AB (r,9 A ,9 B ) = Z, 



2 Pi 

bare *i 



1 + fll ( sin 2 A + sin 2 Q b) + ( sin 2 Q A sin 2 Q B + 1 g m 4 Q A gin 4 + Q(K 6 rl) 

o 64 o 



(7) 

Both expressions agree at order (rero) 2 , and the difference at order (wo) 4 is the monopole-hexadecapole contribution 
which has not been included in (|^) . This comparison illustrates the perturbative nature of the potential derived in 
p3[ , and the fact that the corresponding multipolar contributions are implicitly re-summed in expression (p|) . 



III. GENERALIZATION TO A POLYION OF ARBITRARY SHAPE 



The method used in the previous section may be generalized to find the far field electrostatic potential of an 
arbitrary polyion with vanishing internal volume and bare charge Z — / polyion c(s) cPs. In general, the convolution 
solution (||) admits the large distance behaviour: 



r 



r)~Zl B f(r, K ) + O —5- for «r > 1. (8) 



In this expression, the anisotropic part of the potential again factorizes from the r-dependence and is given by 

/(?,«)= / ^ exp(-Kf-B) d\ (9) 

where f denotes a unit vector in the direction of the position r where the potential is computed, and a(s) is the 
surface charge density at point s on the colloid. As expected, the isotropic bare Coulomb potential is recovered in the 
limit of low electrolyte density: 

lim /(?,«) = 1, (10) 

K, — »0 

with Z the bare charge. 

For a circle of radius ro and uniform line charge we find 

f(0,Kr ) = l (nr sm9), (11) 

where 9 is the angle between r and the normal to the plane containing the circle. Finally, for a uniformly charged rod 
of vanishing radius, and length 21: 

smh(nl cos 9) , x 

mm) = — — jj-A 12 

kI cos 9 

where 9 is now the angle between r and the rod direction. 
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IV. TOWARDS CHARGE RENORMALIZATION FOR PLATELETS 



The bare charge of Laponite platelets can be considered to be of the order of Zb ar c = 700 to 1000 negative elementary 
charges 23 124 . At the level of PB approximation, this corresponds to a high electrostatic coupling (0 larger than 



unity |L3 



), where the linearization procedure underlying Eq. (Q) fails. However, a few Debye lengths away from 



the colloid, the potential has sufficiently decreased so that ([]]) is recovered and the one-body potential takes the form 



e 



4>{r,d) = Z cS f(9,nr 

i 

The effective charge Z c s is (in absolute value) smaller than the bare one |34|-|36f| , and within PB theory, saturates to a 
value independent of Zb arc if the latter quantity is large enough. The contribution embodied in /(#, ktq) results from 
the anisotropy of the colloid, and is a priori an unknown function. In the limit of small electrostatic coupling (small 
^barc) where Z ff/Zb a r C — » 1 by definition, the results of section |l| show that f(9,X) = 2Ii (X sin 9) / (X sin 9) . On the 
other hand, at higher Zb arc corresponding to the saturation regime, the functional form of f(9,X) is the signature of 
the effective charge distribution on the platelet and related to possible differences in counter-ion condensation around 
the centre of the discs or in the vicinity of the edges. 

It has been shown within PB theory that in the colloidal limit ktq > 1, highly charged colloids behave as far as 
their far field is concerned as constant potential particles whose value is close to AkT/e (i.e. cj> = 4), irrespective of 
shape (planar, cylindrical, spherical) p9| |. This prescription may be applied to the present case, with the restriction 
that to our knowledge, the linearized Poisson-Boltzmann problem cannot be solved analytically with the boundary 
condition of constant surface potential. However, the potential associated with the constant surface charge boundary 
condition provides a reasonable first approximation providing the correct qualitative variation of the effective charge 
with physico-chemical parameters. The corresponding anisotropic contribution to the large distance field is that 
computed in section O and the calculation of Z e g follows from the knowledge of 4>{r = 0). From the analytical 
expression reported m [[ll| , we obtain the saturation value of Z e g 

^saturation _ rr _ r 2nrQ 

Z cS — Z sat — -J--. - f r ■ (,-LoJ 

t B 1 - exp(-Kr ) 

It is noteworthy that for a typical salt concentration corresponding to ktq = 1, Z sat ~ 100, which is an order of 
magnitude smaller than the bare charge. This asymmetry guarantees that the effective charge is in the saturation 
regimes where it no longer depends on the value of Z^ arc |34|,|39| (the precise experimental determination of Zbare 
is therefore not required). This saturation picture is expected to be reliable in a 1:1 electrolyte (PB mean-field 
approximation generally fails in electrolytes of higher valence), when there is a clear scale se par ation between the 
Bjcrrum length and the size of the charged object (this constraint being fulfilled for Laponite) |4Ct |. 

The above prediction may be checked by inserting Z sat into the analytical expression giving the LPB force between 
two parallel platelets Jll| (configuration depicted in Fig. ||): 



r F z = iirl B 2 
kT r sa ' 



f Jl{x) ~ exp j-(Vro)y^ 2 + « 2 rg j dx (14) 



and comparing the results with the Monte Carlo simulations of Meyer et al |22]| , who considered the situation of 
vanishing salt. In this limit, where the quality of our prescription for Z e g is expected to deteriorate, we get from 
( |l3| ) Z sat — » 2r /iB — 42. The comparison is displayed in Fig. || which shows a good agreement. Neglect of charge 
renormalization effects lead to an overestimation of the force by more than two orders of magnitude [more precisely, 
by a factor (^b a rc/^s a t) 2 i see the difference between dashed and dotted curves in Fig|^|, which points to the prime 
importance of such a phenomenon. 



V. TENTATIVE INVESTIGATION INTO THE PHASE BEHAVIOUR 



At this point, it is interesting to investigate at least qualitatively the difference in the phase behaviour between 
spherical and disc-like charge stabilized colloids. For spheres, the fluid-solid transition driven by repulsive electrostatic 
interactions favors the isotropic fluid upon addition of salt in the solution Q| . The opposite is observed for Laponite 
dispersions, where an increase of the ionic strength lowers the density where the "solid" phase appears [E 12 1. Given 
that van der Waals interactions are believed to be irrelevant in the corresponding parameter range [12JJ this effect 
(illustrated in Fig. 0), appears at first to contradict standard DLVO phenomenology, where a screening of electrostatic 
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repulsion (decrease of Debye length is expected to promote the simple fluid phase (sol), as happens for spheres. 
However, the effective charge of an arbitrary charged object is generically a growing function of salt concentration 
[see for instance expression (|l3|), which increases with k], as a result of an enhanced screening of colloid/micro- ions 
attraction which diminishes the amount of counter-ion "condensation" . This increase of the effective charge favours 
the solid phase and is thus antagonistic to the abovementioned decrease of the range of electrostatic repulsion driven 
by the decrease of k^ 1 . The competition between the increase of the amplitude of repulsion and the decrease of its 
range upon adding salt is a possible scenario to interpret the phase behaviour of discs, within the DLVO picture and 
without any attractive interactions. More quantitative results are given below. 

In the following analysis, we treat spheres and discs on equal footings (both are considered to have the same radius 
tq). From the above analysis, we consider the following pair potential for discs: 



7 2 „ A Ii(/wosin0yt) li(Kr sm9 B ) e" 



V AB (r, e A ,e B ) = zt at l B 4 ±v -. ~' " ". , (15) 

Kro sin &a kvq sin Ob t 



where Z sa _ t is given by Eq. (|13[); its counterpart for charged spheres is the standard isotropic DLVO expression [KG 32 
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VAB(r) = Z L i B ^—J — (16) 

where the saturation value of the effective charge has been derived analytically in |3^] as a function of salt concentration 
[we again consider colloids with a high bare charge such that Z r coincides with its saturation value]. 

Given that the phase behaviour of particles interacting through Yukawa-like potentials has been extensively explored 
by computer simulations |f42]-f44|| , we can readily obtain the melting density corresponding to a given salt concentration 
for the potential (fl£|). The results (corresponding to spherical colloids) are shown in Fig. |^. A related procedure may 
be used to estimate qualitatively the position of the melting or freezing line for discs, by assuming that at low clay 
density (typically where the fluid/solid transition takes place for Laponite suspensions), the rotational motion of the 
platelets occurs on a shorter time scale than the translational one, so that these objects feel an effective potential 
resulting from the angular average of the expression given in Eq. (|15|), namely 

VlT aS >) = (VAB{r,6 A ,6 B )) 0Af)B = - ( sm0 A d0 A f sm6 B d6 B V AB (r,6 A ,6 B ) 

4 Jo Jo 



,. . , . 2 — * 

72 



cosh(ftro) — 1 \ e~ 

T 

o 



Z L£b4[ ^ ) — • (17) 



In doing so, we obtain an isotropic Yukawa potential where the energy scale -the term in parenthesis in (fl7|)- reflects 
the original anisotropy of the pair potential. In the standard DLVO potential ( |i"(f ) for spheres, the energy scale 
[the term exp(«xo)/(l + Kro)] also depends on Debye length, but has a different physical origin and results from the 
exclusion of micro-ions from the interior of the spheres. Making use of the numerical Yukawa phase diagram [^2[-p[, 
we obtain the melting line represented in Fig. || for the averaged potential (|l^), where the threshold density decreases 
when the salt concentration (or equivalently k) increases, at least for nd < 4. At the same level of description, 
spherical colloids show the opposite behaviour, for all values of nd (we also emphasize that for the spherical colloids 
used in the experiments of Monovoukas and Gast Jtl[, t he equivalent of the melting line reported in Fig. || is in 
excellent agreement with its experimental counterpart |p9| ). 

For both spheres and discs, the prefactor of the Yukawa term exp(— nr)/r is an increasing function of k [see Eqs. 
( |l6| ) and JT^)]. The subtle interplay between this increase and the decrease of Debye length is nevertheless able to 
produce a "re-entrant" melting line for discs only. In the limit of low nr , both expressions ( |l6| ) and ( |i"7| ) become 
£ B Z^ at exp(— nr)/r; the melting lines for spheres and discs however do not merge in this limit in Fig. ^, due to the 
difference in the limiting values of Z sat for both geometries. 

This approach consisting in averaging the two-body platelet potential over angular degrees of freedom predicts 
a qualitative change in the phase behaviour of platelet systems for a reduced density p* = pd 3 of order 1, which 
corresponds to a clay mass fraction of 8% for Laponite with diameter d — 300 A , see the upper x- label of Fig. |. 
This density is approximately 4 times higher than the maximum density delimiting the fluid and solid regions in the 
phase diagram of Laponite suspensions (Fig. || and references [[|[l2| ) , so that the present approach does not allow a 
quantitative comparison with experiments. It is however noteworthy that a density p* ~ 1 is much smaller than the 
isotropic/nematic coexistence density for uncharged plates (p* ~ 4, [EEEq]). At this stage, it is impossible to be more 
specific concerning the nature of the "solid" phase supplementing the fluid one at high densities. This question will 
be addressed in the following section by computer simulations. 
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VI. MESOSTRUCTURE: MONTE CARLO SIMULATIONS 



We have implemented standard Monte Carlo simulations with typically N = 500 platelets interacting through the 
potential (|l5j) where Z sat is given by Eq. (|l3|). This allows to test the validity of the approach proposed in section [v|, 
where the initial anisotropic potential was mapped onto the isotropic Yukawa function (17). A typical run consisted 
in 10 5 cycles (both random displacement and rotation of the N particles). 

For different values of salt concentration and platelet density two diagnostics were used to characterize the simulated 
samples. First, the centre-to-centre pair distribution function g(r) was computed. Second, the orientational ordering 
was quantified through the statistical average of (twice) the second Legendre polynomial P(i/j) — Scos 2 ^ — 1 at a 
given centre to centre distance r, where ip is the angle formed by the normals to the two discs. The orientational pair 
correlation function denoted g or (r) follows from averaging over all pairs of platelets. 

The corresponding information on the mesostructure is displayed in Figs. ^ and [?] for different densities and two 
values of k. A striking observation is that a slight increase of the density induces an important increase of the structure 
(see the differences between the curves at p* = 0.82 and p* = 0.85 in Fig. [| the same is seen at nd = 3 in Fig. |?j 
when the reduced density changes from 1.07 to 1.10). Such values of p* lie close to the threshold p* = 1 estimated in 
section |y| The associated orientational distributions (not shown) are structure-less for the densities considered here. 
These results indicate a qualitative change in the structure of the fluid phase upon increasing the density, but it is 
still difficult to characterize the new "phase" emerging without a thermodynamical study. 

We have tested the relevance of the scenario propounded at the beginning of section [v| (competition between an 
increase of the amplitude and a decrease of the range of the pair potential when n -or equivalently n sa it _ increases). 
For a given density (p* — 1.1), the pair correlation functions are monitored for various values of k. The corresponding 
charges change with k, according to Eq. ([l3|). Figure || shows an interesting feature. As expected, the structure is 
most pronounced for the lowest k, and decreases when k increases (see the difference between the curves for nd = 1 
and nd — 2). However, when salt concentration is further raised to nd = 3, the maximum of g(r) and the whole 
structure are enhanced, before decreasing again when nd = 4. The effect evidenced in Fig. ^ is reminiscent of the 
"anomalous" slope of the phase diagram reproduced in Fig. ||, and in qualitative agreement with the simplified phase 
diagram drawn in Fig. |5|. 

Finally, we have performed exploratory runs at higher densities where a nematic phase would be observed in the 
uncharged system. The results of Fig. [)] for a high density p* — 5.0 show a rich local structure both for the g(r) 
and the orientational g r{r), very different from that observed at lower densities. The reference uncharged system 
(q* = 0) with its strong nematic plateau is displayed for comparison. Inclusion of electrostatic interactions preserves 
the long range nematic order, but decreases its strength (the height of the plateau). Such a nematic order was clearly 
absent at lower densities (see Fig. 0). The peaks of g(r) appear correlated to those of g r{i"), which points to the 
existence of oriented micro-domains with parallel platelets, and a higher nematic order than the mean one. The 
complex behaviour of the distribution functions of Fig. bl calls for more thorough investigations at high densities. In 
particular the validity of the large distance expansion (fLa) may be questionable at high p* . 



VII. CONCLUSION 



In this article we devise an electrostatic pair potential for both highly anisotropic and highly charged objects (of 
vanishing internal volume) based on a linearized Poisson-Boltzmann theory. On the example of infinitely thin platelets, 
we show that the anisotropic shape of the object results in an anisotropic Yukawa potential even at large distances. 
The anisotropy is an increasing function of uro and becomes very marked in the colloidal limit ktq > 1. To account 
for the non-linearities in the Poisson-Boltzmann theory, when the bare charge of the platelets is large, we introduce 
an effective or renormalized charge. This approximation yields an effective pair potential for highly charged platelets. 

The addition of salt generically results in two antagonists effects, irrespective of colloid geometry: the range of the 
effective potential decreases, but its amplitude increases. This interplay appears quite subtle, and able to discriminate 
between spherical and disc-like colloids by producing a re-entrant melting line for discs. This phenomenon appears at 
first incompatible with naive DLVO expectation. Moreover, the corresponding threshold density for discs is smaller 
than the isotropic/nematic coexistence for uncharged plates, but still higher than its experimental "Fluid/Solid" 
counterpart. 

Preliminary results of Monte-Carlo simulations for particles interacting through an anisotropic potential provide 
informations on the mesoscopic structure of the dense phase. In particular, we show that the resulting structure is 
very sensitive upon increasing the density. We have also observed a "re-entrant" melting curve. These effects are in 
qualitative agreement with the experimental phase diagram of Laponite suspensions found in the literature, even if 
this topic is still under debate. 
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Our approach suffers from several weaknesses, a) the non-linearities in the Poisson-Boltzmann theory are accounted 
for at the level of charge renormalization, b) in this framework, we expect the effective charge to be better approximated 
with an anzatz where the particle effectively behaves as a constant electrostatic potential object p^j, whereas we 
considered the platelets to be of constant surface charge, c) we work at the level of pair potential and correlation 
effects beyond the mean-field Poisson-Boltzmann approach might play a role, d) the precise role of van der Waals 
forces has not been assessed, and remains obscure for clay platelets. Further theoretical and numerical work is needed, 
but we hope that the arguments presented here provide a first hint into the full problem. At least, our results strongly 
suggest that the combined effect of anisotropy and charge condensation have a significant and non-trivial qualitative 
influence on the phase diagram of highly charged colloids. These features could well be of prime importance in our 
understanding of the thermodynamics of clay suspensions. 

Acknowledgments. We thank J. -P. Hansen, H. Lekkerkerker, P. Levitz, A. Delville, D. Bonn, T. Nicolai, B. Jancovici 
and F. van Wijland for fruitful discussions. 



APPENDIX A: 

In this appendix, we show that the solution of Eq. (jl]) is not given by convolution product (^) for polyions having a 
non zero internal volume. For simplicity, we consider a charged hard sphere of radius ro (which consequently excludes 
micro- ions from its interior), with a uniform surface charge a. The solution of Eq. (Q) takes the well known DLVO 
form 



(r) =Zh 



I + nr Q r 



where Ze = Anr^a denotes the total bare charge. Alternatively, the convolution route of Eq. pi) yields 



0convol(O = h 



Zt } 



a exp(— k\i 
e |r — : 
sinh(ftro) e~ 



s|) 



d 2 s 



nr 



for r > tq. 



(Al) 

(A2) 
(A3) 



The origin of this discrepancy comes from the fact that upon summing up the elementary surface contributions in 
([A2]), the micro-ions are allowed to enter the "interior" region r < r<j. The field created is thus that of the shell with 
charge Z (the colloid) plus that of the "plasma" inside. The corresponding interior charge Z ln is easily computed 
from the expression of the electrostatic potential for r < r , which follows from a simple permutation of r and r in 
expression ( A3|) 



0convolM = Z£ B 



sinh(/tr) e 



ro 



We obtain 



Z\ n — 



-Ze- Kro 



6^(r)d 3 r 



cosh(Kro) 



sinh(fiTo) 



nr 



Gathering results, a straightforward calculation allows to rewrite Eq. (A3) in the form 

sinh(«ro) e~ Kr 



i(r) 



Z£ f 



Kr 



(Z + Z in ) h 



1 + Kr 



(A4) 

(A5) 
(A6) 

(A7) 



so that the DLVO structure of expression (Al) is recovered. 

The argument given here bears some similarities with the proof of the equivalence between two charged hard spheres 
models with a uniform background 47 (note that the equivalent of the background, the micro- ion density, is not 
uniform in the present situation) . Finally, a related remark is that appl icat ion of expression (^]) for sp here s does not 
give the DLVO pair potential ( |l6| ) even if the correct one body potential ( |A i[ ) is used [use of expression (A2) is equally 
incorrect] . A simple way to recover ([jl]) is through the stress tensor (||] . 
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FIG. 1. Definition of the coordinates used in the two body problem 
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FIG. 2. Geometry used to compute the force between two discs in section 
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FIG. 3. Comparison of the effective force between two parallel and coaxial platelets obtained by (N, V, T) Monte Carlo 
simulations in with the general expression for the LPB force supplemented with the effective charge given by Jl3l). h is the 
distance between the two plates (see Fig. fy. The same results are shown neglecting charge renormalization, which amounts to 
considering Z e Q — Zbare (upper dashed curve) . 
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FIG. 4. Schematic phase diagram of Laponite suspensions, reproduced from Ref. MJ, where the fluid-solid transition was 
referred to as a sol-gel transition. The salt concentration on the y-axis is expressed in mol (dm) -3 , and the scale for Laponite 
concentrations on the x-axis is in g.l -1 . The solid under consideration is isotropic, but at higher clay densities, a nematic solid 
may be formed The corresponding diagram for spherical colloids would exhibit a Fluid/Solid line of opposite (positive) 
slope 
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FIG. 5. Tentative phase diagram of highly charged colloids in a 1:1 electrolyte, as obtained for particles interacting through 
a Yukawa potential in p^-^. The equivalent of the fluid-solid transition for spheres is shown for discs [i.e. the underlying 
Yukawa potential is given by ([li)) for spheres and by jl?! ) for discs]. Discs and spheres have the same radius ro = 150 A 
(reasonable for Laponite clays), and d = 2ro denotes the diameter. The bottom a;-label corresponds the dimensionless number 
density, while the upper i-scale converts this quantity in terms of the ratio of clay mass over solvent mass for Laponite platelets 
(this scale is consequently irrelevant for spherical colloids). In both cases, the line shown is the melting curve delimiting a solid 
at high density and a fluid for more dilute suspensions. 
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FIG. 8. Pair distribution function for p" — 1.1. The reduced charges defined as q" = Zsat y/^B/d have been computed from 
( fi"3| ) and are given for completeness for every k. Also shown is the g(r) for an uncharged platelet system at the same density 
(only excluded volume effects). 
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